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ABSTRACT 

The problem of dynamical heating of galactic discs by spiral density waves is dis- 
cussed using the shearing sheet model. The secular evolution of the disc is described 
quantitatively by a diffusion equation for the distribution function of stars in the space 
spanned by integrals of motion of the stars, in particular the radial action integral and 
an integral related to the angular momentum. Specifically, disc heating by a succession 
of transient, 'swing amplified' density waves is studied. It is shown that such density 
waves lead predominantly to diffusion of stars in radial action space. The stochastical 
changes of angular momenta of the stars and the corresponding stochastic changes of 
the guiding centre radii of the stellar orbits induced by this process are much smaller. 
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1 INTRODUCTION 

The shearing sheet (Goldreich & Lynden-Bell 1965, Julian 
& Toomre 1966) model has been developed as a tool to study 
the dynamics of galactic discs and is particularly well suited 
to describe theoretically the dynamical mechanisms respon- 
sible for the formation of spiral arms. For the sake of sim- 
plicity the model describes only the dynamics of a patch 
of a galactic disc. It is assumed to be infinitesimally thin 
and its radial size is assumed to be much smaller than the 
disc. Polar coordinates can be therefore rectified to pseudo- 
cartesian coordinates and the velocity field of the differential 
rotation of the disc can be approximated by a linear shear 
flow. These simplifications allow an analytical treatment of 
the problem, which helps to clarify the underlying physical 
processes operating in the disc. 

In two previous papers (Fuchs 2001a, b, referred here- 
after to as papers I and II) I considered the stellardynamical 
model of a shearing sheet and discussed the unbounded sheet 
and then the dynamical consequences, when inner bound- 
ary conditions are applied. The aim was to give a consistent 
theoretical description of pure swing amplification (Toomre 
1981) as well as exponentially growing modes in the frame- 
work of the same model. 

It is well known from numerous studies that, if a star 
disc is perturbed by a succession of spiral density waves, the 
stars are scattered randomly by the spiral arms and their 
velocity dispersion grows steadily so that the background 
states of the discs are evolving (Julian 1967, Carlberg & Sell- 
wood 1985, Binney & Lacey 1988, Jenkins & Binney 1990). 
Dynamical disc heating has also been demonstrated in nu- 



merical simulations of the dynamical evolution of star discs 
such as by Sellwood & Carlberg (1984) or Toomre (1990). 

In the present paper I discuss disc heating of the shear- 
ing sheet. I follow in particular the theory of difi'usion of 
stars in the two-dimensional action integral space due to 
wave - star scattering developed by Dekker (1976). In sec- 
tion (2) I briefly describe the formal derivation of the difi^u- 
sion equation in the framework of quasi-linear theory (Hall 
& Sturrock 1967) and in section (3) I calculate diffusion co- 
efficients for scattering of stars by swing amplified density 
waves. This kind of spiral density waves has been shown 
to be relevant for disc heating by Toomre (1990), wheras 
mode-like density waves, on the other hand, do not heat ef- 
fectively (Barbanis & Woltjer 1967, Lynden-Bell & Kalnajs 
1972). 



2 DERIVATION OF THE DIFFUSION 
EQUATION 

When describing the disc heating effects of a succession of 
transient spiral density waves, one has to distinguish be- 
tween the long time scales, on which the overall distribution 
function of stars in phase space is evolving, i.e. the dynami- 
cal heating time scale, from the shorter time scales, on which 
the individual density waves develop. As is well known from 
plasma physics (Hall & Sturrock 1967) this concept allows to 
derive from the Boltzmann equation in quasi-linear approx- 
imation a Fokker-Planck equation for the long term evolu- 
tion of the distribution function. I follow here in particular 
the adaption of the formalism to the dynamics of stellar discs 
by Dekker (1976). The phase space distribution function can 
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be expressed with the aid of the action and angle variables, 
Ji, J2 and wi,W2, respectively, 



/( Jl, J2,Wl,W2; t) . 



(1) 



The evolution of the distribution function with time is de- 
termined by the coUisionless Boltzmann equation. 



0, 



(2) 



where H is the Hamiltonian of the stellar orbits, and the 
square bracket denotes the usual Poisson bracket. Following 
Dekker (1976) I define a suitable mean of the distribution 
function by averaging it over the angle variables. 
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dwi / dw2f{Ji,J2,wi,W2;t). 
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(3) 



This averaged distribution function will evolve slowly on the 
long time scale. For the rapidly fluctuating part of the dis- 
tribution function 5f — f— {/) one obtains from the general 
Boltzmann equation (2) 

^ + ^ + [{f),Ho] + [Sf, Ho] + [{/) , S<i>] + [Sf, 5$] = ,(4) 

where the Hamiltonian has been split up as // = Hq + 5^ 
with (53> denoting the fiuctuations of the gravitational po- 
tential. The first Poisson bracket in equation (4) vanishes, 
because {/) depends only on the action integrals. The time 
derivative of the averaged distribution function is expected 
to be much smaller than that of the fluctuating part. Thus 
equation (4) can be cast into a Boltzmann equation for Sf, 
and, if the quadratic term [5f, S^] is neglected, into a lin- 
earized Boltzmann equation. 



[Sf,Ho] = -[{f),5'i>] 



(5) 



dsi 

dt 

which has been widely used to study the dynamics of galactic 
discs. {/) has taken the role of the axisymmetric, stationary 
background distribution. Its slow evolution with time is ne- 
glected on short time scales, but taking the average of the 
Boltzmann equation (2) leads to 



dt 



(6) 



The second and third terms of equation (6) vanish because 
{5f) = (5$) = by definition, so that equation (6) takes 
the form 



+ =0, 



(7) 



dt 

which descibes the long-term evolution of the averaged dis- 
tribution function in quasi-linear approximation. Equation 
(7) is also valid, if instead of action and angle variables other 
variables are used, as long as Ji, J2 are integrals of motion 
and wi, W2 their canonical conjugates. 

The potential perturbation 5$ can be Fourier trans- 
formed with respect to the angle variables. 



5$ = / dh / dl25<^i{3;t)e 



and similarly 

5f=( dh I d/2<5/i(J;t)e'''^™^+'""'"' 



(8) 



(9) 



On the other hand, the Ihs of equation (5) represents the 
total derivative of 5f along stellar orbits in the axisymmetric 
part of the gravitational potential. Thus 



5f 



dt' [{f){3,,-t'),5^[3,,,w,,-t')\ 
j dt'Y,^\ I dh I dl2S^,{J,,;t') 

J to n ' J J 



(10) 



where the integration is to be taken along 'unperturbed' or- 
bits. The indices of the action and angle variables indicate 
that the variables, which are the independent variables of 
the distribution function and the gravitational potential, re- 
spectively, must be chosen according the 'unperturbed' orbit 
starting at Jtp, Wtg and terminating at Jt, Wt. In the next 
section I will apply equation (10) to a succession of uncorre- 
lated swing amplification events of short duration. The typi- 
cal integration interval t — to will be then much smaller than 
the time scale, on which the averaged distribution function 
{/) is evolving and equation (10) can be simplified to 



5f = 



d{f) 

^ dJn 



to 



dt' / dh / dl2S^i{Jt,t') 



Comparison of equations (9) and (11) shows that 



(11) 



^^i^n f dt'SMJut') 

.g»('l(™i,t' ~™l.t) + '2(m2,t' -™2,t)) 



(12) 



Upon inserting expressions (8) and (9) into equation 
(7) it is straightforward to evaluate the Poisson bracket and 
carry out the averaging with respect to the angle variables. 
After some algebra one obtains 



d{f) 
dt 



dh 



dhY.^n^iSfiS'i';) . 



(13) 



where use of the fact has been made that 5<I>_i = 5^1 so that 
5$ is real. Since the potential perturbations are supposed to 
be a succession of short lived, uncorrelated fluctuations, it is 
customary to include into the averaging process an ensemble 
average over these fluctuations. Inserting finally equation 
(12) into (13) leads to a diffusion equation in action integral 
space. 



d{f) 



iv^ d 



dt 2 ^ dJn 

m,n 



d{f) 



dJn ' 

with diffusion coefficients 



(14) 



A,,n = 2 j dh j dhkJu j dt'{^{t)^x{t')) 

.g»('l{™l_t' -1"l,t) + i2{™2,t' -1"2,t)) 



(15) 



Equations (14) and (15) are valid for any gravitational per- 
turbations of the stellar disc with moderate amplitudes 
and short correlation time scales. Very similar relations 
have been derived by Binney & Lacey (1988) in a different 
way under more general assumptions. However in Dekker's 
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(1976) approach the duahty of equations (5) and (7) is par- 
ticularly instructive. 



3 DISC HEATING BY SWING AMPLIFIED 
SPIRAL DENSITY WAVES 

I apply the formalism of the previous section to calculate 
the diffusion coefficients for wave - star scattering by shear- 
ing, swing amplified spiral density waves. The dynamics of 
the density waves are modelled by a shearing sheet made 
of stars. This describes a patch of a thin galactic disc. Its 
centre orbits the galactic centre at galactocentric radius vq 
with an angular velocity Qq- Pseudo-cartesian coordinates 
are defined with respect to the centre of the patch, 



r-ro,y = ro{9 - Qot) , 



(16) 



where r and denote galactic polar coordinates, respec- 
tively. As explained in paper (I) the equations of motion 
of the stars are derived from the Hamiltonian 



Ho = ir' + ^rlie - Qof - 2Ano{r - ., 
or alternatively 

nr r 1 ^ r2 ^ n2 2 

Ho = kJi + — Ja - -iZoro , 



(17) 



(18) 



where A and B denote Oort's constants. The resulting orbits 
are simple epicyclic motions, which can be written as 



_ + y_sm^i,2; = ™2- 



2B 



cosiiJi , (19) 



with K = \/—4:^loB the epicyclic frequency. Ji is the radial 
action integral of an orbit. J2 denotes the integral J2 = 
y + 2flox of an epicyclic orbit and is related to the angular 
momentum of a star as J2 = {r^9 — roQo)/ro- As can be seen 
from equation (19) the guiding centre radius of the orbit is 
given by 



J2 

-25 



(20) 



wi — Kt and u)2 ~ %J2t are variables canonical conjugate 
to Ji and J2, respectively. The radial and circumferential 
velocities are given by 



2B 



u = \/2kJi coswi, V = \J2kJ\ sin w\ . 



(21) 



where v is defined relative to mean shearing velocity y = 
—2Ax. The dynamics of a shearing sheet made of stars has 
been studied extensively by Julian & Toomre (1966) (cf. also 
Toomre 1981) and is discussed at length in paper (I) using 
strictly Eulerian coordinates. The principal result is that 
the wave crests of the density waves, which appear in the 
disc, swing around following the mean shearing motion of 
the stars. While the waves swing around their amplitudes 
are amplified transitorily and then die away. In paper (II) 
it is shown that also mode - like, quasi-stationary solutions 
can be constructed for the shearing sheet by introducing an 
inner reflecting boundary in the disc. Such kind of density 
waves is of no concern in the present context, however, since 
they hardly heat the disc at all. 



The perturbations of the gravitational potential of the 
shearing sheet are customarily Fourier analyzed as 



(5$ : 



dkx 



k^ „ „ -j- k^ 



-2B 



sm wi 



\/2kJi 

+ kyW2 — ky — — COS wi 

2B 



(22) 



The functional dependence on the variable W2 is of the form 
as in equation (8) and I use in the following the wave number 
I2 instead of ky. The dependence on the angle variable wi 
can be easily adapted to the form of equation (8) by taking 
an inverse Fourier transform of equation (22) with respect to 
wi. The resulting Fourier coefficients are inserted then into 
equation (15). There is, however, a difference of the mean- 
ing of the variable W2 in this section from the meaning of 
the corresponding variable in the previous section. W2 mea- 
sures now the drift of the guiding centre of the orbit in the 
y-direction and is thus not an angle variable. Averaging the 
distribution function with respect to W2 is impractical and 
I consider in the following the evolution of the distribution 
function in ( Ji , J2)-space, i.e. the distribution function in- 
tegrated over wi and W2 , respectively, which is conceptually 
slightly different from the average (3). This leads to expres- 
sions for the diffusion coefficients of the form 




dl2 I df'eH-'i"'i+*^ 



. r , , I , / / 2,7i . , , ^2kJi 



(23) 



Note that the diffusion tensor is symmetric by construction. 

In paper (I) it is illustrated how the shearing sheet re- 
sponds to internal and external perturbations by developing 
density waves with shearing crests and amplitudes ampli- 
fied while the waves swing around. Fig. 1 shows the distri- 
bution of amplitudes of the potential perturbations as func- 
tion of wave numbers k^,ky, respectively, calculated using 
the formulae of paper (I) (section 11.4) for the case when the 
swing amplification mechanism is fed by white noise. Toomre 
(1990) discusses the source of the noise and argues convinc- 
ingly that swing amplified white noise explains exactly the 
behaviour of the star discs in his numerical simulations or 
that by Sellwood & Carlberg (1984). The distribution of am- 
plitudes, which represents the superposition of many short 
lived shearing density waves, is quasi-stationary and can be 
modelled for positive wave numbers ky empirically as 



-k^Or (fcy-fcyO)" 



5<I>k = l>e 



(24) 



and continued at negative wave numbers ky as (J^-k = 'J'I'k. 
The parameters are estimated as fcxo = 1-5 fccrit, fcyo = 
0.5 fccrit, (Tk^ — 0.7 fccrit, aud (Tky =0.1 fccrit for the case 
A = —B = ^f^o. The critical wave number is defined as 
fccrit = K^/(27rGEd) with G the constant of gravity and Ed 
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Figure 1. Response of the swing amplifier to continous white 
noise. The quasi-stationary distribution of the amplitudes of the 
potential perturbations of the star disc, ^l^^p, is shown as func- 
tion of radial and circumferential wave numbers fcx 
negative wave numbers ky the distribution of amp 
tinues as (5<E>_k = S^i^. Each spike represents a superpostion of 
swing amplified density waves all travelling at a constant wave 
number ky along the fcx abscissae at a speed of fcx.eff = 2Afcy. 
The wave numbers are given in units of fccriti whereas the ampli- 
tudes are normalized arbitrarily. The adopted parameters of the 
disc model are A = —B = if2o and Q = 1.4. See paper (I) for 
further details. 



the surface density of the disc, respectively. Using this para- 
metric model the quadratures in equation (23) can be carried 
out explicitely. Each spike in Fig. 1 represents the super- 
postion of swing amplified density waves all travelling at a 
constant wave number ky along the fcx abscissae at a speed 
of fcx,cff ~ 2Aky with an amplitude which can be modelled 
as 



^ 



(25) 



In equation (25) fcJf denotes the original radial wave num- 
bers of the waves and t'" the time, when the waves were 
launched, ^^.m i^i^^"^) parameterizes the excitation rate of 
the waves and the peak amplification factor of the swing am- 
plifier (cf. paper I) . Since white noise is a stationary random 
process, <E>j;in k2(*'") evenly distributed over wave number 
space and t'". Each of the density waves is only correlated 
with itself. This implies, if the autocorrelation function in 
equation (23) is considered, restrictions of the wave numbers 
and the time intervals to 

k]^ + 2Al2{t' -t'") - ^ 0, and 

k^ + 2 Ah (t - t'") - fcx = so that 

fci - fcx - 2 Ah (t - i') = . (26) 

This constraint and equation (25) lead to an autocorrelation 
function in the form of 

{5$k,,i2(t>'l>ki,i2W) 

= |$o|'<5(fcx-fcx-24;2(t-t')) 

(fex-fexp)^ (fcx-fcxp)^ ('2-''yo)^ 

•e ^"'L e '"kx e "^y , (27) 

where |$o|^ is a normalization constant. The time integral 
of the autocorrelation function according to equation (23) is 



given by [h > 0) 

dt'5{k'^ - fcx - 2Ah{t - t'))c'['i"+'=^4-'2](t'-t) 



_1_ i[;i^+,2 4j2]^x^ 



2Ah 



(28) 



I consider next in equation (23) the integration with respect 
to wave number fcx. 



(fcx-fcxo)'' 

27^ 

dfcxC e 



2vr(jt e 



/2.7i . ii-t 



(29) 



and similarly 



(fcx-'=xn)-' 



dk'^e 

oo 



— v27ro"kxe 



/ 2Ji . / , il« 



In equations (29) and (30) terms of the kind 



2 /2Ji 
'^kx 



(30) 



(31) 



are neglected as quadratically small. This is justified, be- 
cause the majority of the stars have epicycle sizes smaller 
than the critical wave length Acrit, the typical spacing be- 
tween spiral arms (Julian & Toomre 1966, cf. also paper 
I). The epicycle size is determined by ^2J\I whereas 
(Jkx oc fccrit = 27r/Acrit (cf. Fig. 1). Similarly terms of the 
kind 



2 



f— v 

will be neglected, because 



OC — < 1 . 



2^^2 fcx cff ^orb Torh 



(32) 



(33) 



where Tacc denotes the duration of disc heating by a single 
spike of the potential fiuctuations, which is only effective 
close to the peak in Fig. 1, whereas Torb is the orbital pe- 
riod of the stars. In this aspect the disc heating mechanism 
described here is rather impulsive. Next I consider the inte- 
gration in equation (23) with respect to l\, 



(34) 



where the upper row refers to Dii, the middle row to D\2 = 
D21, and the lower row to -D22, respectively. The results are 
given by delta functions and derivatives thereof. 
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{—5"{wi — w'l) 
iS'iwi - w[) 
S{wi — w'l) 



(35) 



The next step is the integration with respect to the angle 
variable wi, which leads for the diffusion coefficient D22 to 



by multiplying the diffusion coefficients by a factor of two. 
Assembling all results leads a diffusion tensor of the form. 



)Jl 



(42) 




dwiS{wi — 'u;i)e 



{sin TiJi —sill ti;^ ) -j^ 



/ Np-''2 'ili ^ (costui-cos w[) 



(36) 



■ I J / /^ ~-il2 OP (cos wi —cos wi ) 

I / dwid [wi — Wi)e ^ 2B V ij 



-I / 2-/1 , . . / . 

^0 v — ;t— (sm KJi — sin J 



2B 

and for Dn to 



—- — sm -u;! — fcxo \ / cos 

IB V ^ 



(37) 



dwiS" (wi — w'i)e 



zfexQ Y/ — j7— (sin — sill ti;^ J 



/ \ — 1(2 — 2S — ^"^^^ "^1 — cos w^) 



'0 



., , ., /2Ji 

th cos Wi - Jfcx„ A / sm Wi 

Id V ^ 



/ . y^ini . , /2Ji 
I — — sm w-i + ifcxQ y cos -u;! 




(38) 



The integration over w'l gives simply a factor of 27r for Dii, 
whereas the integration over the trigonometric functions in 
equation (37) leads to the result 



D12 = D21 = . 



Integrating equation (38) with respect to w'^ gives 



/2 2kJi 2 



2J^\ 



4B2 

The final step is the integration with respect to I2, 

('2-fcyO)^ 
71 



(39) 



(40) 



t2 



1 



7J- (^;? j_ j,2 2J]_ 

1 



(41) 



where the upper rows refer to the diffusion coefficient Dn 
and the lower rows to D22, respectively. There is a formal 
divergence at ^2 = in the second term of Dn on the Ihs of 
equation (41). I have chosen to ignore this and have replaced 
the integral by a saddle point approximation. The reason is 
that the model of disc heating used here (cf. equation 25) 
becomes unphysical for small circumferential wave numbers 
I2, because the density waves approach then the WKB limit 
and become long lived, so that they do not heat the disc 
effectively. Due to the symmetry of the distribution of am- 
plitudes (24) with respect to k the effect of density waves 
with negative wave numbers ky can be taken into account 



with Do = 87r2 |<I>opo-^ (Jk. 



fcaoM. 

The diffusion equation takes the form 



d{f) 
dt 



(43) 



where the overhead tilde means the the Ji -dependence of 
Dii has been written separately. For this kind of disc heat- 
ing by transient spiral density waves no correlation between 
the diffsion in radial action and angular momentum space 
is found. The diffusion equation (43) is highly non-linear, 
because the diffusion coefficients depend on the distribu- 
tion function (/) themselves. In particular, the effectivity 
of swing amplification of spiral density waves depends crit- 
ically on the value of the Toomre stability parameter Q — 
K(Ju/(3.36GEd), where Uu denotes the radial velocity disper- 
sion of the stars. Thus, when the disc heats up, the ampli- 
tudes $0 of the density waves and the diffusion coefficients 
(42) will decrease and the disc heating rate slows down to 
zero. Numerical simulations of the dynamical evolution of 
galactic discs (Sellwood & Carlberg 1984, Fuchs & v. Lin- 
den 1998) have shown that this can happen on compara- 
tively short time scales, if the discs are left uncooled. Only 
if the discs are cooled dynamically by adding stars on low 
peculiar velocity orbits, the spiral density wave activity can 
be maintained at a constant level despite the rising velocity 
dispersions. In that case Do ~ const, and a simple solution 
of the diffusion equation (43) is found by separation of vari- 
ables in the form 

(/> = — ' 



ici + iq^t)Jc2 + 



exp - 



Ji 



+ 



ci + 



C2 + ^t 



(44) 



with arbitrary constants Ci and C2. 



4 DISCUSSION AND CONCLUSIONS 

The scattering of stars by shearing, short lived density waves 
leads to independent diffusion of stars in radial action - an- 
gular momentum space. This diffusion process has various 
implications. The radial action integral is related to the pe- 
culiar velocities of the stars as 



" +4B5" 



Thus a distribution function of the form 



exp-Ji/ ( ci + 



(45) 



(46) 



implies a exp —u^ /2a^{t) dependence of the velocity distri- 
bution with a predicted radial velocity dispersion of (Tu(t) = 
^yK.{cl + bllt/2). Such a rise of the velocity dispersion with 
time fits ideally to the actual age-velocity dispersion relation 
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observed in the solar neighbourhood (cf. Fuchs et al. 2001 
for a recent review of the observational data) . Unfortunately 
the diffusion coefficients cannot be estimated quantitatively, 
because the constant I'l'oP, which parameterizes the white 
noise, is not known a priori. But judging from the shape 
of the heating law (46) wave - star scattering of the kind 
discussed here might well have played an important role in 
the Milky Way disc. Whether this is the only disc heating 
mechanism is still a matter of debate (Fuchs et al. 2001). 

The diffusion of the guiding centre radii of the stellar 
orbits can be estimated from the ratio of the diffusion coef- 
ficients, 



D22 



Jl k2 



(47) 



As shown above KDwtj (2 Ji) is proportional to the square of 
the radial velocity dispersion, a\{t). According to epicyclic 
theory (cf. equation 20) 



(48) 



so that the dispersion of the guiding centre radii of the stellar 
orbits is given by 




t = 



-2B 



,(49) 



+ 



Using the parameter values estimated above equation (49) 
implies ^ (xg) — 0.2(Ju/f2o, if a flat rotation curve is as- 
sumed. A star with an age of 5 Gyrs like the Sun has typi- 
cally in the sol ar ne ighbourhood a velocity dispersion of 50 
km/s. Thus (xf) = 400 pc — 0.05 ro in the solar neigh- 
bourhood, if a local angular velocity of fio ~ 26 km/s/kpc 
is adopted. This confirms the conclusion of Binney & Lacey 
(1988) that the diffusion of guiding centre radii driven by a 
rapid succession of spiral density waves is rather small. This 
assumes that disc heating by transient spiral density waves 
is the only disc heating mechanism. However, if the Sun 
has indeed drifted from its birth place nearly 2 kpc radially 
outwards to its present galactocentric radius as suggested 
by Wielen, Fuchs, & Dettbarn (1996) and Wielen & Wilson 
(1997), this would mean that there must be other dynam- 
ical heating mechanisms of the galactic disc. It was shown 
by Fuchs, Dettbarn, & Wielen (1994) theoretically and by 
numerical simulations that Spitzer-Schwarzschild diffusion 
of stars due to gravitational encounters with massive molec- 
ular clouds, for instance, leads to a much more pronounced 
diffusion of the guiding centre radii of the stellar orbits, even 
if the mechanism does not heat the disc effectively. 
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